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ABSTRACT 


This study represents performance tests of the third-order turbulence 
closure for predictions of separating and recirculating flows in 
backward-facing step. Computations of the momentum and temperature fields in 
the flow domain being considered entail the solution of the time-averaged 
transport equations containing the second-order turbulent fluctuating 
products. The triple products, which are responsible for the diffusive 
transport of the second-order products, attain greater significance in 
separating and reattaching flows. Formulations have been made for those 
products mentioned above by developing each corresponding transport equation. 

A low Reynolds number model associated with the viscous effect in near-wall 
regions is developed and incorporated. 

The computations are compared with several algebraic models and with the 
experimental data. The predication has been improved considerably, 
particularly in the separated shear layer. 

Computations are further made for the temperature-velocity double 
products and triple products. Although there is not any appropriate 
experimental data to compare, the present computations are compared to the 
results obtained by solving the several existing algebraic correlations. The 
agreement between these sets is shown to be quite reasonable. 

Finally, several advantages have been observed in the usage of the 
transport equations for the evaluation of the turbulence triple products; one 
of the most important features is that the transport model can always take the 
effects of convection and diffusion into account in strong convective shear 
flows such as reattaching separated layers while conventional algebraic models 
cannot account for these effects in the evaluation of turbulence variables. 


NOMENCLATURE 


2 

wall skin friction coefficient, t. ,/(l/2)pU Tlll 

W IN 

constant used in turbulence model 

coefficients for the pressure-heat flux term in the transport 
equation of <u^0> 

coefficients in the algebraic expressions for <u.e> 
coefficients for the pressure-stress and the near-wall 
low-Reynolds number model of the transport equations of 

<u i u jV 

coefficient for the dissipation rate in the transport 
equations of <u i u j u | c > 

coefficient for the determination of turbulent viscosity 
(= 0.09) 

coefficient for the pressure-heat flux term in the transport 
equation of <u.u.0> 

• J 

coefficient for the dissipation rate in the transport 
equation of <u.u.e> 

* J 

specific heat of fluid 

wall static pressure coefficient, (P w - Pj N )/(l/2)pUj N 

depth of the channel (Y + H) 

diffusion rate of Reynolds stresses, <u..u.> 

diffusion rate of triple-velocity products, <u^UjU^> 

diffusion rate of the heat flux, <u^0> 

diffusion rate of the triple products of velocity and 

temperature, <u_.Uj0> 

step height 


t ' 


P ijk,l 


P ijk,2 


P ije,l 


ije,2 


turbulent kinetic energy 
pressure fluctuation 
mean pressure 

production rate of Reynolds stresses, <u.u.> 

* J 

production rate of triple velocity products, <u.u.u.>, due to 
mean strain rates J 

production rate of triple velocity products, <u.u.u.> , due to 

1 J K 

interaction of the Reynolds stresses with their gradients 
production rate of <u.e> 

production rate of <u.u.0> due to mean strain rates 

* J 

and temperature gradients 

production rate of <u.u.e> due to Reynolds stresses 

* J 

and the heat fluxes interacting with their gradients 

molecular Prandtl number 

turbulent Prandtl number 

wall heat flux per unit area 

mean temperature 

inlet stream temperature 

fluctuating velocity in x-direction 

mean velocity in x-direction 

inlet stream velocity 

fluctuating velocity in y-direction 

mean velocity in y-direction 

reattaching point 

normalized coordinate = (x - x )/x. 

r r 

Cartesian coordinates 
normal distance from the wall 

height of the channel upstream of the step 
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Greek Svmbols 


a 

thermal diffusivity 

4 ij 

Kronecker delta 

c 

dissipation rate of turbulent kinetic energy 

c ij 

dissipation rate of Reynolds stresses, <u ^ u j u | c > 

c ijk 

dissipation rate of triple-velocity products, <u.u.u.> 

i J Is 

c ie 

dissipation rate of the turbulent heat flux, <u..e> 

c ije 

dissipation rate of third moments of velocity and 
temperature, <u.u.e> 

* J 

e 

fluctuating component of temperature 

V 

dynamic viscosity 

v t 

2 

turbulent viscosity (= C^pk /e) 

V 

kinematic viscosity 

P 

density of the fluid 

T 

shear stress 

d) t • 

pressure-strain redistribution for Reynolds stresses, 
< Vj > 

*1 j .w 

near-wall effects in the pressure-strain redistribution for 
Reynolds stresses, <u.u.> 

' J 

^ijk 

pressure-stress correlation for triple-velocity products, 

<u.u.u, > 
i J k 

*ijk,w 

pressure-stress correlation of the near-wall low-Reynolds 
number effect for triple-velocity products, <u..UjU k > 

‘he 

pressure-heat flux correlation for <u..e> 

4> . 

ie,w 

near-wall effects in the pressure-heat flux correlation for 
<u i e> 

d>. 

MJ6 

pressure-heat flux correlation for the triple products, 

<u u e> 
i J 
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Subscripts 


i ,j ,k,ft.,m,n 

w 

IN 

Symbols 
< > 


tensor notations 
wall values 

inlet station (x/H = -4) condition 


time-averaged values 
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INTRODUCTION 


Separation and recirculation of flow is encountered in a vast array of 
engineering applications. This phenomenon of flow separation and recircula- 
tion associates itself with higher turbulence levels which not only render the 
flow greater analytical complexity but also result in augmenting its heat 
transfer and momentum aspects considerably. 

Current research in this area is mainly limited to experimental measure- 
ments of the hydrodynamic characteristics along with the heat transfer 
coefficients in the flow domain. Computation of the temperature field in the 
flow domain entails the solution of the time averaged temperature equation 
containing the turbulence fluctuating velocity-temperature products. These 
second moments or double products are correlated in terms of Boussinesq's 
isotropic viscosity and the turbulent Prandtl number in order to close and 
solve the temperature equation. However, as pointed out by Launder and 
Samaraweera , 1 the second-moment model is more preferable because turbulent 
interactions which generate the turbulence stresses and heat fluxes can be 
treated exactly. They applied an algebraic second-moment turbulence closure 
to heat and mass transport in thin shear flows and obtained good results. 

Algebraic turbulence scalar-velocity models containing the curvature and 

2 

gravitational effects have been proposed and used by Tahry et al. and 
3 

Gibson for predicting the heat fluxes in turbulent shear flows and boundary 
layer flows. 

The third moments, which are responsible for the diffusive transport of 
the second moments, attain greater significance in recirculating flows, as 

4 

pointed out by Chandrsuda and Bradshaw ; these also have not been given much 
attention and, therefore, need to be examined in greater detail in order to 
fully understand and comprehend the turbulent transport mechanisms. 
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Owing to the relative complexity commonly encountered in recirculating 
flows, experimental measurements of the second and third moments, hence, pose 

5 

a formidable challenge. Ota and Kon report uncertainties of up to 50% in 
their measurements of second moments in the reattachment region using conven- 
tional hot wire instruments; similar measurements were made by Seki et al. 6 
7 8 

Vogel and Eaton and Adams et al. made extensive measurements of the 

thermal and mass fluxes using LDAs. Experimental information on the behavior 

of the triple products of scalar-velocity fluctuations is almost nonexistent. 

g 

Andreopoulos and Bradshaw , however, present the variation of the triple 
products for a flow over a boundary layer using constant temperature 
anemometers and Antonia 1 ® presents those for a plane jet. 

The study undertaken and presented here solves the flow into a backward- 
facing step geometry. After obtaining a good agreement for the hydrodynamic 
parameters, the thermal variables are solved. The transport equations of the 
turbulence scalar-velocity products, <u..e> and <u i Uj6>, are 
developed and solved to obtain their profiles in recirculating and 
redeveloping flow regions beyond the step. 

MATHEMATICAL MODELS 


The transport equation for the Reynolds stresses is given as 


sir (“k < u i“j>) ■ p 1 j - c ij + ♦jj + *ij,w + °1J 


( 1 ) 


where P.., <(>.. , and 4>... , respectively, represent the production, 

1 J 1 j 1 J ,w 

pressure-strain, and the wall correlation for pressure-strain rates. These 
are all defined by Launder et al. 11 

The diffusion rate, D.., contains the terms given as follows 

' J 
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( 2 ) 


<u.p> <u.p> 

D-. = - zr~ [<u.u.u,> + — 4.,+ - J — 
ij ax k 1 i j k p jk p ik 


- V 


3<U-U.> 3U. 3U. 

(— **— + <u. T^> + <U. — ->)] 

' ax k j ax. i axj n 


Since the diffusion rate of the Reynolds stresses is governed by the gradient 
of the third moments <u i u j u | C > » it becomes necessary to evaluate the 
third moments accurately. 

The third moments, <u ^ u j u | c > t can be evaluated by formulating their 
transport equation as 


au k au, au* 


ax7 < u i <u i u jV> = " (<u i u jV iT + <u j\V ST + <u k u i V ajJ) 


(i) 


a <u kV 


+ (<u i u j > sr~ * <u j v + <u k u i > -it~> 


a<u^Uj> 


(II) 


au k u i 3u i u i 3u i u « 

- (<u.u. ■ „ — ■ > + <U.U. — -> + <U, u. - " J > ) 

' U ax t J k ax 4 k i ax^ • 

(III) 

■ p (<u i u j ^ - <o j u k ^ * <u k u i 

( IV) 

a 3u k 3u « 

+ <u i u j ax^ [v ( ax t + ax k )>] 


( 3 ) 


a 3u i 3u i 

+ < Yki^ [v 


au. au 

+ <u, u. - 2 - [v ( — ^ + — -)>] 
k i ax £ 1 v au t axj ; J 

(V) 
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In order to account for the viscous effects that are predominant near the 
wall, a low-Reynolds number modification needs to be incorporated in the 
transport equations of the third moments. 

The theory behind this low-Reynolds number modeling is based upon the 
existence of the viscous sublayer, the buffer layer and the fully turbulent 
core in the near-wall region. The profiles of the energy dissipation rate c 
are accounted for in all three layers, and the pressure-stress correlation is 
modified accordingly. 

Finally, the transport equation for < u i UjU k > containing the low- 
Reynolds number effect is given as 


<u iU .u k >) « P. jki+ P ijk2+ 4> ijk ♦ * ijkfW 


+ c ijk + D ijk 


( 4 ) 


where 


P i j k 1 = production due to mean strains 

3U k 3U. au. 

■ -C g* <u 1 u j u l > + <VkV 5 * <u k u 1 V 


(5) 


P ij k ,2 = Production due to Reynolds stresses interacting with their 
gradients 


3<u.u.> 3<u,u.> 3<u.u •> 

= " {<U kV ~ 3X + <U ^ U o > + Jt ) 


11 9X„ 


J a ax. 


( 6 ) 


j k + ^ijk w = Pressure-stress term including the low-Reynolds 
number effect. 
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(7) 


<u i u iV k 3/2 ak 1/2 2 

- S ~k tc + max t c Y w hr • 2v <rh- } ]} 

X 


= Dissipation due to viscous action 
- <Vr ! (4 1j 4 4 jk + 4 ki> * k ’ /? 


and = Diffusion rate 


“ (v 5^ <u i u jV) 


(8) 


(9) 


The values of the constants are listed in Table 1. 

1 2 

The model given by Eq. (7) was originally tested by Amano et al. 
where it was shown that the addition of the bracketed term remarkably improves 
the prediction of <u.u^Uj c > in the near-wall region. 

Table 1. Recommended values for the constants used 
in turbulence modeling. 


C 9 

c 

y 

c 

YW 

C 

CY 

c ie,l 

c ie,2 

C IT 

C 2T 

C e Y 

C 0c 

1.0 

3.0 

8.0 

0.10 

3.2 

0.5 

0.313 

0.156 

6.0 

0.10 


The Reynolds-stress closure is consolidated by incorporating the equa- 
tions of k and c and the transport equations of the third moments, <u.u.u.>, 

1 J K 

into the Reynolds-stress equations. 

The temperature field is obtained by solving the time-averaged tempera- 
ture equation given as 



rt_ £L] 

p rj ax, J 


(10) 
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where Pr t = turbulent Prandtl number (= 0.9). 

In order to obtain the second moments of temperature, <u^e>, their 
transport equation is formulated and various terms in it are then approximated 
through closures. 

The transport equation for <u.e> after neglecting insignificant terms 
can be written as 


zr~ (U. <u.e>) = - [<u.u.> + <u.e> — -1 

ax. v j i ' L i j ax. u j ax. J 

J J J 

(I) 

+ afr + v > 577 <u i 6> - <u i u j 6> ] 

J J 

(II) 


au. 


P 30 

+ < E zr~> 
p ax. 

(ill) 


ae 8u i 

(a + v) <5 *j 
(IV) 


( 11 ) 


which can be written as 


iJT < yV» 


= P . + D . „ + <b . - c. 

10 10 M0 K 


where 


( 12 ) 


P. fl = Production rate of <u.0> 

D.* = Diffusion rate of <u.e> 

iw i 

4 >. 0 = Pressure-heat flux effects 
and c. Q = Dissipation of <u.0.> 

I" 1 

P i0 ( Term (I)) needs no further approximation since it is explicit in 
its character. 
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The diffusive rates which are comprised of the triple products 
<u.u.e> are modeled by decomposing the third moments into the second 

’ J 

1 3 

moments as shown by Launder, 


v d<u.e> 

<u.u j e> = - 0.15 (2 *) <u A > —^] 


Thus, the diffusion is now given as 


a<u.e> 3<u.e> 

°ie - iJT [<- ♦ v) -sj- ♦ < Uj V 


The pressure-heat flux term can be expressed as 


♦le ' ' c ie,i k <u i 6> + c ie,2 <u » 9> 5^ 


The dissipation rate, c._, is assumed to be negligible in accor- 

1 v 

dance with the assumptions made in the Reynolds-stress closures. Launder and 

Samaraweera 1 also propose a near-wall correction term to be 

1 9 f w 

added to the pressure-heat flux term, <t>. Q , which is given as 

au. au 3/2 

<t>. ^ = [- o.io <u.e> - 0.02 <u.e> {4 — 1 - — ^}] - — m 

ie,w L k 1 i ax a ax^ J c x n v 


where x^ is the distance from the wall, 
n 

Algebraic expressions to model the second moments are also developed by 
1 4 

Launder which after neglecting the gravitational effects can be expressed 


k aT k 9U i 

<u.e> = C 1T - <u. u.> rj— + C, T - <u u e> — 

i 1 T c k i ax k 2T c k ax k 
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In order to evaluate the third moments of temperature, <u.u.e>, it 

* J 

is necessary to formulate their transport equation. This is given as follows: 


^ < u k <u i u j 9> > - 


au. au, 

' [<u l“jV 5^ + <u j“k 6> ^ + <u k u i 6> a^ 

(I) 


- [< Vj 3^ <V ‘ <u k e>)> + <u i e 3^ ( "j“k - <u jV )> 


* <u j e 5^ 'Vi - <u k u i >)>] 

(in 


- [<^ 2E-> + <!!j! 5 e_>] 
L p ax.. P ax i J 

(in) 


au. 


- [<U • U • (a “ — )> + <U .0 ( 

1 i j ax k v ax k / i ax k K P ax k 


au. 


)> + <u .6 t~" ( M T- 1 

• -I * A ^ V 


j ax. 'p ax. 


(IV) 


(18) 


where 

Term (I) = Production due to mean strain rate and temperature gradients 
Term (II) = Production rate due to the interaction of the Reynolds 

stresses and their gradients with the heat flux components 
<u..e> and <Uj©> . 

Term (III) = Pressure-heat flux effects 
Term (IV) = Diffusion and dissipative effects due to molecular viscosity. 

In closing the above equation. Term (I) needs no further approximation 


since it is rather explicit in its form. Term (II) can be rearranged and 



3<u.e> a<u.u.> a< u. ii . > 

Term (II) - < Vj > — - + <u,e> -J— * < Uj 6> — — - 


( 19 ) 


The quadruple term is assumed to be Gaussian and can be split up as 


<u.UjUj < e> = <u.Uj>»<u^e> + <UjU^>*<u .e> + <u^u.>*<Uj©> 


( 20 ) 


Differentiating Eq. (20) with respect to x k and substituting the result in 
Eq. (19) yields 


Term (II) = 


a<u.e> a<u.e> 

- <u,u,.> -r: <u, u .> — 


jk ax. 


k i ax. 


a<u.u.> 

- <m e> — 1 
k ax,. 


( 21 ) 


The pressure-heat flux effects of Term (III) can be approximated by first 
forming the Poisson equation for fluctuating pressure and then approximating 
the term by a simplified closure. 

The Poisson equation for fluctuating pressure is 


l an a 2 2 au au 

~p i^7 ' " J ( Vm - <u l u m > ^ + ~ax7 ST } dx t 

JL si m Si m 


( 22 ) 


After multiplying u p e and time averaging with suitable approximation, 
this may be given as 


U 0 

2E_> = 

p 


au 

- C l <u n u l e> k + C 2 <u m u r> 9> JT 

m 


(23) 


where C-j and are arbitrary constants. 

Based upon this approximation. Term (III) is thus correlated as 


u.e 


u .6 


Term 


(III) = - <— — - <-*•— = c. <u.u.e> r 

v ' p ax. P ax. l i j k 


au, au, 

+ c 2 [<u j u k e> aT + <u k u i e> ax7^ = c e Y <u i u j e> k 


( 24 ) 
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where the second term in Eq. (24) is merged into Term (I) and the coefficient 
C adjusted accordingly. (The value of has to be determined 
by parametric testing.) 

Term (IV) contains the diffusive and dissipative effects due to molecular 
viscosity. This may be rearranged as 


Term (IV) = 


<“ + 2 »> ^ < Vj e>1 


- a 


3U,U. 63 U.U. 

[2< |f- + < V- >) 

3x k 3x k 3X Z 


2 * 

3U. 3U .6 U.3 U.6 

- v [2<— y-> + < ^->) 


3X k 3X k 


3X, 


2 * 

3U • 3U.6 3 u.e 

■ v [2< jt + <u i — r >] 

3x k 3x k J 3X.. Z 


(25) 


The first term in Eq. (25) represents the laminar diffusion while the 
rest of the terms express the dissipative effects. The terms with the 
asterisk (*) consist of second derivatives of the second moments and are 
assumed to be negligible. The effect of neglecting the terms is adjusted by 
manipulating the value of the empirical coefficients used in the modeling. 

The dissipative effects of Term (IV) in Eq. (25) are further assumed to be 


3c 


r . . = c - (<u,e> — ) 
lje 6c c ' k 3x . ' 


(26) 


The final form of the transport equation for <u.u.e> is now written as 

* J 


^ (U k <u.u.e>) = P ije l + P ij0 2 - D. . e - 4>. . e - c. j0 


(27) 
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where 


t au. au. 

p ije,i * - [< VjV * <u j u k e> + ‘Vi 9 * ^ 


a<u.e> 


a<u.e> 


o = - [<U^U,,> + <U L U.> -z^— + <^6> J ] 


ije,2 


j k ax. 


a< u.u.> 


k i ax k ' '“k" ax k 


D 1je ■ «“ + 2 »> ^ ‘Vj 9 * 5 


d> . = - r r <u . u . 6> 

v ije ©y k i 3 


and 


: . . = C f <u, e> 

13© ©c k k ax. 


In addition to the transport equation for <u.u.©>, the existing 

* %) 

algebraic correlations are: 


1 . Launder 


13 


<u.u.e> 

* s) 


2k 


a<u.©> 


a<u^©> 


• - 0-10 r (< 1 , V -a * <u jV IT") 


2. Wyngaard and Cote 


15 


< u 1 “j e > - - ° ' 5 f < <u iV -ii-> 


a<u.e> 

c i 


3. Donaldson et al. 


16 


<u 


iUj e> = - 0.10 (— 


..2 a<u.e> a<u.©> 
4k__ , t . j , 


ax. 


SOLUTION PROCEDURE 


(28) 

(29) 

(30) 

(31) 

(32) 


(33) 


(34) 


(35) 


The computation was performed using a 62 x 62 variable grid mesh with the 
grid expanding linearly at the rate of 2% in the axial direction and at 3% in 
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the transverse direction from the step wall. The details of this adoption is 
discussed in the next section. 

Computations of the transport equations described in the preceding 

1 2 

section is achieved by using a finite volume method. The iterative 
procedure is terminated when the maximum value of the relative residual 
sources of U, V, and mass balance falls below 1%. However, the computations 
of the triple products are terminated when the relative residual sources fall 
below 3 x 10' 8 for <u.u.u.> and 5.0 x 10 -9 for <u.u.0>. 

This complete process of solving the momentum, temperature and their 
related turbulence products equations takes about 60 minutes of CPU time on a 
UNIVAC 1100 computer. For the solution of the transport equation of the 
temperature and temperature associated moments, boundary conditions had to be 
specified at the solid wall and the outflow section of the channel. 

For the temperature equation, a constant heat flux of 130 W/m was applied 
along the wall downstream of the step. This heat flux was introduced into the 
solution domain by supplementing the source term at the wall adjacent cell for 
the temperature equation with this heat flux. 

The wall boundary condition for <ue> and <ve> is based upon the fact 
that at the wall, 

<u6> = <V0> = 0.0. (36) 

Therefore, the value of <v0> at the wall adjacent node is fixed by the 

interpolation between the wall and the node next to the wall adjacent node. 

The value of <ue> at the wall adjacent node is set equal to minus twice the 

13 

value of <v©> at this node, as pointed out by Launder. 

Since the scalar-velocity products <u.u.0> fall to zero at the 

* J 

wall, their near-wall values, hence, have been accordingly made very small 
such that they tend to zero at the wall. 
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12 

As mentioned in Amano et al. the third moments are fairly stable and 
seem to be affected very little by the choice of the wall boundary treatment. 
At the outflow section which is located at a distance of 60H from the step, a 
parabolic boundary condition is applied thereby causing the streamwise 
gradient of all the parameters (except T) to vanish. 

The temperature distribution at the outflow section is prescribed in 
accordance with Kays and Crawford 17 wherein 


dj 

dx 


li 



PC UD r 

p c 


(37) 


RESULTS AND DISCUSSION 

Figure 1 depicts the flow configuration inside a backward-facing step 
geometry. The shear flow separates at the edge of the step and forms a 
separated shear layer which causes a part of the fluid to recirculate and the 
rest to continue to flow downstream contributing to the stabilization of the 

flow through the formation of a boundary layer in the redeveloping region. 

1 8 

Experimental data of Driver and Seegmiller has been used for compari- 
son with the computational results. 

Figures 2 and 3 compare the skin friction coefficient and pressure 
coefficient, respectively along the wall downstream of the step with the 
experimental data. Computations have been performed with several different 
grid systems ranging from 32 x 32 to 62 x 62. It is shown in both figures 
that a fair grid independent status is attained with the grid system finer 
than 52 x 52. Thus the grid of 62 x 62 has been adopted in this study. It is 
further noticed that agreement between the computation and the data is fairly 
reasonable. 


17 - 


Figure 4 compares the mean velocity profiles and Figs. 5, 6, and 7, 

2 2 

respectively, compare the Reynolds stresses <u >, <v > and <uv> with the 
experimental data. The agreement with the experimental data is generally 
satisfactory. 

In Figs. 8-11 triple velocity products are compared with the experimental 

1 8 

data of Driver and Seegmiller. The results are also compared with several 
algebraic models in the literature. As shown in these figures, the levels of 
the triple-velocity products are always higher in the shear layer with 
employing the transport models than with the algebraic correlations. Thus, 
the transport model gives more realistic results in comparison with the data. 
The reason for this is that the convective and diffusive effects are accounted 
for through the reattaching shear layer with the transport model whereas these 
effects are mostly neglected in the algebraic correlations. 

Since the physical properties of the fluid are assumed to be invariant 
with temperature, hydrodynamic aspects of the flow are independent of its 
thermal characteristics. Therefore, the temperature along with its second and 
third moments are obtained using the converged stored values of the hydrody- 
namic parameters. 

Figure 12 shows the nondimensional ized temperature profiles downstream of 

1 9 

the step. These results, comparing the data taken by Eason and Vogel, 
show very satisfactory agreement. 

Figures 13 and 14 show respectively the nondimensional ized heat flux 
<ue> and <ve> profiles downstream of the step. Results obtained by the 
solution of their transport equation (Eq. 12) are compared with the results 
obtained by the algebraic correlation (Eq. 17). It is seen that the overall 
agreement between the results obtained by the transport equations and the 
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algebraic equations is quite reasonable. It is further noticed that the 
agreement of <v©> is very good in the shear layer but closer to the wall the 
algebraic model fails to diminish the levels of <ve>. On the contrary, the 
algebraic models tend to make <v©> go to infinity. The transport equations, 
however, with the proper imposition of the boundary conditions give more 
realistic results. 

Profiles of the triple products of velocity and temperature <u.u.e> 

^ J 

are shown in Figs. 15-17. Results by solving the transport equations for 
<u.u.©> (Eq. 27) are compared with results obtained by the algebraic 

* J 

13 15 

correlations for <u-u.©> by Launder , Wyngaard and Cote, and 

* J 

Donaldson et al., 1 ^: Eqs. (33), (34) and (35), respectively. 

There is no experimental data available for validation and comparison 

owing to the relative complexity involved in the measurements of double and 

triple products (<u.©> and <u.u.©>) inside wall bounded recirculating flows. 

Figure 15 shows the nondimensionalized profiles for <uu©>. Predictions 

by the three algebraic correlations seem to be very close to one another. The 

predictions by the transport equation are the highest particularly in the shear 

layer. Although the three algebraic correlations predict similar results, the 
1 3 

model of Launder predicts the highest values and that of Donaldson et 

al. 16 predicts the lowest. However, the difference among the predictions is 

very small, particularly near the reattachment region and beyond. 

Figure 16 depicts the profiles of <uv©>. The prediction by the 

algebraic models and the transport equations are very close. The model of 

Donaldson et al.^ gives the highest values whereas the model of Wyngaard 
1 5 

and Cote gives the lowest. Predictions by the transport equations are 
somewhat intermediate. 
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Figure 17 shows the profiles for <vve> obtained by the transport 
equation (Eq. 27) and the above mentioned algebraic correlations. The predic- 
tions by the model of Donaldson et a!. 16 predominate over the others, and 

the transport equations give the smallest values. The predictions by the 
1 5 

model of Launder show smaller than those by the model of Donaldson et 
al., 1& but larger than those of Wyngaard and Cote's 15 model. The differ- 
ences between the predictions become smaller with increasing distance down- 
stream from the step. 

It is generally seen that the predictions by the algebraic models and the 
transport equations have similar trends and the difference between their 
predictions becomes smaller in the reattachment and redeveloping regions. 
However, the near-wall values of all the models are rather high. It is 
possible to treat the near-wall values of <u.u.e> through the incorpora- 

* J 

tion of a suitable low-Reynolds number model in their transport equations. 

But it is not easy to modify the algebraic equations to accommodate this. It 
is anticipated that the transport equations of the third moments with a 
low-Reynolds number model would be superior to other algebraic correlations. 

CONCLUSIONS 

The second-and third-order closures for the hydrodynamic characteristics 
are extended to the scalar (temperature) field wherein transport equations for 
the fluctuating velocity-temperature moments are developed and solved for the 
flow over a backward-facing step. 

Prediction for <u^8> and <u.Uj0> obtained by solving their transport 
equations are compared with algebraic models and the superiority and 
potentiality of the transport equations for these double and triple products 
are demonstrated. 
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FIGURE CAPTIONS 


Figure 1 

Figure 2 

Figure 3 
Figure 4 
Figure 5 
Figure 6 
Figure 7 
Figure 8 
Figure 9 
Figure 10 
Figure 11 
Figure 12 
Figure 13 
Figure 14 
Figure 15 
Figure 16 
Figure 17 


Solution domain -- Backward-facing geometry with the separating, 
recirculation and redeveloping regions. 

Skin friction coefficient along the bottom wall downstream on the 
step. 

Pressure coefficient along the bottom wall downstream of the step. 
U-velocity profiles downstream of the step. 

<uu>-prof iles downstream of the step. 

<vv>-profiles downstream of the step. 

<uv>-profiles downstream of the step. 

<uuu>-profiles downstream of the step. 

<uuv>-profiles downstream of the step. 

<uvv>-profiles downstream of the step. 

<vvv>-prof iles downstream of the step. 

T-profiles downstream of the step. 

<ue>-profiles downstream of the step. 

<ve>-prof iles downstream of the step. 

<uu©>-prof iles downstream of the step. 

<uve>-prof iles downstream of the step. 

<vve>-prof iles downstream of the step. 
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Figure 1 Solution domain -- Backward-facing geometry with the separating, 
recirculation and redeveloping regions. 
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Figure 2 Skin friction coefficient along the bottom wall downstream on the 
step. 
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Figure 3 Pressure coefficient along the bottom wall downstream of the step 



Figure 4 U-velocity profiles downstream of the step 
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Figure 5 <uu>-profiles downstream of the step 




EXPERIMENT 
(Driver & Seegmiller) 

COMPUTATION 



Figure 6 <vv>-prof lies downstream of the step. 
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Figure 7 <uv>-prof iles downstream of the step. 




Figure 8 <uuu>-prof lies downstream of the step 



Figure 9 <uuv>-prof lies downstream of the step 



Figure 10 <uvv>-prof iles downstream of the step 




Figure 11 <vvv>-profiles downstream of the step 









TRANSPORT MODEL 

ALGEBRAIC MODEL 
(Launder & Samaraweera) 



Figure 13 <ue>-prof i les downstream of the step 
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ALGEBRAIC MODEL 
(Launder & Samaraweera) 



Figure 14 <ve>-profiles downstream of the step 






TRANSPORT MODEL 

ALGEBRAIC MODELS 
— Wyngaard 

Owen & Launder 

Donalson 



Figure 15 <uuQ>-prof iles downstream of the step 







